Size effect in fracture: roughening of crack surfaces and asymptotic analysis 
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Recently the scaling laws describing the roughness development of fracture surfaces was proposed 
to be related to the macroscopic elastic energy released during crack propagations. On this basis, 
an energy-based asymptotic analysis allows to extend the link to the nominal strength of structures. 
We show that a Family- Vicsek scaling leads to the classical size effect of linear elastic fracture 
mechanics. On the contrary, in the case of an anomalous scaling, there is a smooth transition from 
the case of no size effect, for small structure sizes, to a power law size effect which appears weaker 
than the linear elastic fracture mechanics one, in the case of large sizes. This prediction is confirmed 
by fracture experiments on wood. 

I. INTRODUCTION 

In solid mechanics, an essential scaling problem is the effect of the structure size on its nominal strength. This 
effect is particularly important in the case of quasibrittle materials which are characterized by the existence of a 
large fracture process zone containing many damage microcracks. Materials as different as concretes, mortar and 
rocks, some composites, toughened ceramics andpwood belong to this category. BazantEij has shown that in this case, 
contrary to what happens for Weibull's statistical, the size effect is linked to the very existence of the development of 
such a large microcracked zone, which implies stress redistributions and stored energy release. 

On the other hand, in quasibrittle materials, damage has^a strong influence on the local deviations of the main 
crack through its elastic interactions with the microcracksHia. Consequently, the roughness of the fracture surfaces 
can be considered as an inheritance of the damage process and it is naturally tempting to correlate their morphology 
(and especially their fractal properties) to macroscopic mechanical properties such as fracture energy or fracture 
toughnessErt^l. ,— ■ 

In this paper, on the basis of a link recently established!!^ between the roughening of crack surfaces and the 
macroscopic energy released during crack propagation (this work is summarized in Section Q), we propose an energy- 



based asymptotic analy sis in Section III which allows to extend this link to the nominal strength of large (Sec . [II A) 



and small (Sec. [II B) structures. An approximate size effect valid everywhere is proposed in Sec. Ill C . This 



prediction is shown to be in agreement with experimental results obtained on wood Sec. IV. Finally, we discuss 
material-dependent properties in Sec. |v|. 

II. SCALING LAWS OF CRACK SURFACES AND ENERGY RELEASE RATE 

The statistical characterization of the fractal morphology of fracture surfaces is nowadays a very active field of 
research—It is now well established that-these. surfaces pfjpx, very different types of materials (from ductile aluminium 
alloysoE3 to brittle materials like rocHHEj'ES or woocE3E!l), exhibit self-affine scaling properties in a large range of 
lengthscales (seeO for a more detailed account of experiments). Moreover, in addition to this_self-affine character, 
recent studies focussed on the complete description (3D) of the morphology of crack surfacesE3E3 have shown that 
the scalins-iaws governing the crack developments in longitudinal and transverse directions are different and material 
dependent! 2 !! 

Let us consider the development of a fracture surface from a straight notch of length L with zero roughness. 
The mean plane of the crack surface is defined as (x, y) where the x axis is perpendicular tottlic directiop_of crack 
propagation and the y axis is parallel to this direction. For two quasibrittle materials (granitec3 and woodEj), it has 
been found that the fluctuations Ah of the height on the fracture surfaces, estimated over a window of size / along 
the x axis and at a distance y from the initial notch exhibited anomalous-scaling properties which are quite similar 
to those obtained in some models of nonequilibrium kinetic rougheningc3 ! E3: 

where £(?/) = By x / Z depends on the distance y to the initial notch and characterizes a crossover length along the x 
axis. For length scales smaller than £(y), the surface is self-affine, and characterized by the local roughness exponent 



1 



Q oc . This self-affine character is observed in most experiments and a local roughness exponent Q oc ~ 0.8 is reported in 
all cases. Hence, it has been suggested that this local roughness exponent might be a universal index, i.e. independent 
of the fracture mode and of the materialE-3. 

According to Eq. ([!]), the magnitude of the roughness increases as a function of the distance y until the self-affinc 
correlation length £(y) reaches the system size L. This happens at a certain distance y sa t = (L/B) z from the notch: 
£{y ^ llsat) = L. Thus, the first growth regime of the roughness {i.e. for y <C y sa t) is followed by a stationary regime 
(for y ;g> y sa t) where the magnitude of the roughness remains constant and where the global roughness {i.e. measured 
over the system size L) is driven by the global roughness exponent £: Ah{L, y 3> y sa t) ~ L ■ The main consequence of 
an anomalous scaling [Eq. ([[])] is that, in this stationary regime, the magnitude of the local roughness {i.e. measured on 
windows I -C L) is not only a function of the window size I but also of the system size L: Ah{l, y 3> y sa t) ~ ^ !oc L^~^ loc . 

Experimental results obtained on quasibrittle materials have shown the following values for the global roughness 
exponent: Q — 1.2 for graniteta and 1.35 and 1.60 for wood (pine and spruce respectively )c9. Note that the exponent 
z (called the "dynamic exponent^-} and the prefactors A and B seem also to be material-dependcntnJ. 

The Family- Vicsek scaling lawE-1, where the height fluctuations Ah scale as: 

can be seen as a particular case of the anomalous scaling [Eq. (|l|)] where C = Ooc- In this case [Eq. (|J)] , the magnitudes 
of the local and global roughnesses at saturation, respectively Ah{l,y ^> y sa t) ~ ^ !oc and Ah{L,y ^> y sa t) ~ L^ 1 "", 
are driven by the same roughness exponent (the local roughness exponent Cioc)- Furthermore, the local roughness 
is independent of the system size. This is drastically different from what happens in the case of an anomalous 
scaling. It is worth noticing that both the Family- Vicsek and, the , anomalous scaling laws were used to describe 
the roughness development fti the same granite fracture surfacetZnlS, and that more accurate results were obtained 
assuming anomalous scalingcS. 

The anomalous and Family- Vicsek scalings_were recently shown to be linked to drastically different mechanical 
behaviors in terms of elastic energy releasee!! Within the framework of an equivalent linear elastic problem, a 
fracture criterion, linking the elastic energy release rate G at the macroscale and the fractal nature of the crack at the 
microscale, was proposed: G = 2 r yip(y)/L. In this fracture criterion, 7 is the so-called specific surface energy and the 
ratio ip{y)/L can be considered as a "roughness factor". As a matter of fact, it is precisely the ratio of the length of a 
virtual crack front over its projected length L (the system size, i.e. the specimen width). This virtual front, parallel 
to the initial notch, is assumed to be rough only out of the average fracture plane (x,y), and this roughness can be 
defined from the functions Ah{l,y) at a fixed position y from the initial notch. 

According to an anomalous scalpg [Eq. ([[])] , the estimate of the real length ip of the virtual front (which corresponds 
to the length of a self-affine curveO) leads to the following expressions of the energy release rate: 



G R {Aa)^2 7 {\ lif -I (3) 

£ 2(C ~°° c) if Aa > Aa sat 

where the crack length increments Aa and Aa sa t defined from the initial notch correspond respectively to the crack 
positions y and y sa t defined in Eq. (Q). In the zone where the roughness grows, i.e. for Aa <C Aa sa t, the fracture 
equilibrium leads to an energy release rate function of the crack length increment Aa [Eq. (g)]. The subscript R in 
Gr emphasizes the fact that-the resistance to fracture growth is similar to the behavior described by a resistance 
curve (usually called i?-curveO). Note that the square root terms in Eq. (^) are dimensionless and correspond to the 
roughness factor ip{y)/L. The term l a is the lower cutoff of the fractal range of the virtual front {i.e. the characteristic 
size of the smaller microstructural element relevant for the fracture process). When the crack increment is large, 
i.e. for Aa ^> Aa SQt [Eq. (^)] which corresponds to the saturation state of the roughness, the resistance to fracture 
growth becomes independent of the crack length increment because the self-affine correlation length has reached the 
system size: £(Aa » Aa sat ) ~ L. Introducing the crossover length Lc = {lo 1 ^^ 10 " / A) 1 /^^ 1 "^, the resistance to crack 
growth for large crack length increments becomes: 




G R {Aa » Aa sat ) ~ G RC ~ 2 7 Wl + — (4) 



where the subscript C in G R c emphasizes that the resistance to crack growth has reached an asymptotic or critiml 
value. The main consequence of the link between fracture mechanics and anomalous roughening of fracture surfacesli2l 
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is the size effect on the criticai resistance to crack growth [Eq. (0)]. This size effect is the resuit of the dependence of 
the maximum magnitude of the roughness at saturation on the structure size L. As shown in Fig. [l], the size effect 
affecting the critical resistance Grc exhibits two asymptotic behaviors separated by the crossover length Lc- For 
small structure sizes (i.e. L <C Lc) where the roughness of the fracture surface is weak, there is no size effect, and 
Grc — 27. For large structure sizes (i.e. L 3> Lc), corresponding to an important fracture roughness, the critical 
resistance evolves as a power law : Grc ~ L'»~^ oc . 

In the case of the Family- Vicsek scaling [Eq. (||)] , on the contrary, the link between roughnening of fracture surfaces 
and material fracture propertiestij reduces to a resistance to crack growth independent of the crack increment Aa and 
of the specimen size L: 



Gc(Aa)^2 7 ^l+^-A-) (5) 

Thus, an anomalous scaling accounts for an i2-curve behavior and a size effect on the critical resistance to crack 
growth while a Family- Vicsek scaling reflects a purely elastic brittle fracture behavior. 



III. SIZE EFFECT ON THE NOMINAL STRENGTH 



We now extend the connection just summMizedliij to the nominal strength of structures. 

Within the framework of Bazant's theorycJ, the size effect (for two-dimensional problem) can be described from 
geometrically similar structures of different sizes (with geometrically similar initial cracks or notches) by introducing 
a nominal stress: 
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where P is the external load applied to the structure (considered to be a load independent of the displacement), L is 
the characteristic size of the structure, and d is any length of the structure (for instance, as shown in Fig. || in the 
case of TDCB specimens, d corresponds to the ligament length). When P = P u which corresponds to the ultimate or 
maximum load, <jn is called the nominal strength of the structure. 

On the other hand, at the maximum load P u , the elastic energy release rate G (obtained at a constant load P u 
or cn) must be equal, according to the linear elastic fracture mechanics (LEFM), to the critical resistance to crack 
growth G RC : 
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where the complementary energy W* characterizes the energy stored in the structure. This energy W* can be written 
in the following form: W* = o~n 2 L d 2 f(a)/E where / is a dimensionless function characterizing the geometry of the 
structure and a — a/d is the relative crack length. Thus, when the maximum load is reached, the nominal strength 
of the structure can also be written as: 
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where g(a) = df(a)/da corresponds to the dimensionless energy release rate function. 

Usually, the main problem with Eq. (^|) is to determine the relative crack length a = a/d for which the maximum 
load P u is reached. Let us consider exclusively a structure of "positive geometry", i.e. dg(a)/da > under load 
control. 

If the material exhibits a purely elastic brittle behavior such as the one obtained for a Family- Vicsek scaling [Eq. 
(g)], the relative crack length increment a at the maximum load corresponds to a a = a /d (a being the length of the 
notch or initial crack). Thus, in the case of a Family- Vicsek scaling, replacing the resistance to crack growth Gc by 
its expression in [Eq. ^) and a by a Q in Eq. (|^), we obtain the following size effect relation: 
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where m is a proportionality coefficient between d and L (d = mL). In Eq.(||), the term under the square root is a 
constant, which means that the nominal strength of structures evolves as <7jv ~ L^ 1 ^ 2 . This is in agreement with the 
expected LEFM size effect, i.e. the size effect of a purely elastic brittle fracture behavior. 

In the case of an anomalous scaling, the problem appears more complicated. Indeed, under the conditions evoked 
in the previous case (i.e., structures of positive geometry and under load control), if the material exhibits an i?-curve 
behavior [Eq. (g)], the crack length increment Aa sa t (limit of the i?-curve) is actually the limit of stability. Thus, 
the relative crack length at the maximum load can be defined as a = a + where 9 = Aa sa t / d. Hence, in the case 
of an anomalous scaling, which reflects the fracture behavior of quasibrittle materials, the knowledge of the evolution 
of the crack length increment Aa sa t as a function of the structure size is the key of the size effect problemj—However, 
the dependence between Aa sa t and the specimen size does not appear clearly from the roughness analysisEfrEj. .— ■ 

In order to resolve this problem - which is the central point of this paper -, a possible way suggested by BazantO 
consists in considering that the failure of a quasibrittle material is not only characterized by the specific surface 
energy 27 (related to the actual crack surface), but also by a critical damage energy release rate Gd per unit volume 
of damaged material {i.e. per unit volume of fracture process zone). Thus, one can assume that failure at the maximum 
load is obtained for the energy balance: 

G d V FPZ = 2 7 A r (Aa sat ) (10) 

where the volume of the fracture process zone can be estimated as Vfpz = L Aa sa t 2 /n with LAa sa t the projected 
crack surface and Aa sat /n the height of the process zone (where n is assumed constant, i.e. independent of the 
size L). The surface A r (Aa sa t) corresponds to the real fracture surface produced during the crack advance Aa sat . 
Note that the fracture criterion [Eq. dig)] actually corresponds to an equivalent linear elastic problemE3c3 where the 
effective size of the fracture process zone at the maximum load is assumed equal to the crack length increment for 
which the resistance to crack growth does not follow the ii-curve [Eq. (||)] but remains constant and equal to the 
critical resistance: G R {Aa sat ) = G RC [Eq. (§)]. 



In the following, on the basis of the fracture criterion defined in Eq. (10), the asymptotic values of Aa sa t are 
estimated for large and small structures in order to obtain the nominal strength respectively for large and small-size 
asymptotic expansions. 



A. Large-size asymptotic expansion of the size effect 

As previously mentioned, the square root terms in Eq. (|^) and Eq. (^) correspond to the ratio of the virtual 
crack front length ifj(Aa) over its projected length L. Moreover, from Eq. (Q), relative to the roughness satu- 
ration regime, one can obtain the maximum length of the virtual crack front since ?p(Aa 3> Aa sat ) ~ ipmax — 

Thus, for large structure sizes, i.e. L> Lc, the asymptotic value of the real crack surface 
A r produced for a crack advance Aa sa t can be estimated from ip m ax as: A r ~ (3 L Aa sa t (L/Lc)^ - ^ 1 "" where j3 is a 
constant (function of the scaling exponents) . Substituting A r in Eq. ([To]) yields the expression of the crack length 
increment: 

Aa sat = c*P i—j (11) 

where c* = n2^/Gd is a material-dependent length scale. Thus, for large structure sizes, the relative size of the 
fracture process zone (i.e. 9 = Aa sat /d) is expected to evolve as a power law 9 ~ L^~^ lo ^.L, Moreover, from the 
values of the scaling exponents obtained in the roughness analysis of quasibrittle materialst2lc2l, it appears that the 
relative size of the process zone becomes negligible when the system size increases: lim 9 = for L —> +00. In other 
terms, in large structures, the process zone is expected to lie within only an infinitesimal volume fraction of the body 
and so lim a — a for L — » +00. Note that this result is in agreement with Bazant's assumptiono. Hence, the 
dimensionless energy release rate function g(a) being generally a smooth function, we may expand it into a Taylor 
series around a = a Q and Eq.(||) thus yields : 
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are all constants. Equation ( |13[ ) provides a large-size asymptotic series expansion of the size effect because the terms 
containing non zero powers of L in denominator vanish for L —> 00. Note that Eq. is expected to diverge for 
structure sizes L — > as shown in Fig. ||. A first-order asymptotic approximation at large sizes can be obtained by 
truncating the series after the linear term: 
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(15) 



The main consequence of the anomalous roughening in the case of large structure sizes is that the nominal strength is 
expected tp-decrease as on ~ £-V 2 +(f-0<><!)/ 2 ^ ag s i 10wn m pig ||) This re sult disagrees with the size effect proposed 
by Bazanto where the nominal strength of large structures decreases as on ~ L~ x / 2 which is the theoretical size 
effect of LEFM. The difference originates in the fact that, for an anomalous roughening, the critical resistance to 
crack growth Grc [Eq. (Q)] is expected to evolve as a power law Grc ~ L^^ 1 "" for large structure sizes while in 
LEFM, the critical resistance Grc is assumed to be constant {i.e. independant of the specimen size). Hence, the size 
effect on the nominal strength obtained for an anomalous roughening is weaker than the size effect in LEFM. 



B. Small-size asymptotic expansion of the size effect 



In Bazant's theoryEj, no size effect is expected for small structure sizes (L — » 0); this is the domain of the strength 
theory. A possible justification is that, in small structures, the fracture process zone fills the whole volume of the 
structure and hence, there is no stress concentration and, as a consequence, failure occurs with no crack propagation. 

Such an argument can be also obtained from the link between anomalous roughening of crack surfaces and material 
fracture properties. Indeed, in small structure sizes {i.e. L <C Lc), the roughness being negligible, the virtual 
crack front length tends to its projected length L. This implies that, for small structures, the actual crack surfaces 
produced during a crack advance Aa sat are not so different from the projected one: A r ~ LAa sa i. Hence, substituting 
A r ~ LAa sat into the fracture criterion [Eq. (|lC|)] yields the crack length increment: Aa sa f = n2j/Gd = c* for 
L <C Lc- In other terms, the effective size of the process zone tends to the material length for small structure sizes. 
Thus, when the material length c* = d — a a (Fig. U), the fracture process zone occupies the entire ligament of the 
structure. _ 

On the basis of Bazant's theorjc3 and in order to obtain a small-size asymptotic expansion of the size effect, let us 
now introduce a new variable and a new function: 
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(16) 



The function (p(a , rf) corresponds to the dimensionless energy release rate function of the inverse relative size of the 
process zone. Substituting Eq. ( |l^ ) into Eq. (||) and expanding (p{a ,rf) into Taylor series around the point (a o ,0) 
since lim r\ = when d or L — > 0, yields the nominal strength: 
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(17) 



(18) 



where (pi{a o ,0) = d(p(a o ,0)/dr], (p 2 {a o ,0) = d 2 (p{a , 0)/cV, and c 2 = ip 2 {a , 0)(p(a o , 0) 2 /{2ipi{a a , 0) 2 ), c 3 
^ 3 (a o ,0)(;9(a o ,0) 3 /(6v3i(a o ,0) 3 ), and, 
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cap = 



/ 2-fE 
ip(a O7 0)c* 



Ln=C* 



mcpi(a Ol 0) 



(19) 



are all constants. Equation (|18|) provides a small-size asymptotic series expansion of the size effect and is plotted 



in Fig r 

theory 



When L — > 0, the nominal strength tends to a constant (i.e. cap as expected in the case of a strength 
1)" but diverges from the asymptotic behavior of the size effect obtained in the case of large structure sizes [Eq. 



C. Approximate size effect 

Now, the main problem consists in interpolating between the large-size [Eq. (|l3|)1 and the small-size [Eq. (list)] 
asymptotic aeries expansion in order to obtain an approximate size effect valid everywhere. The theory of intermediate 
asymptoticsE3 is not easily applicable in our case. Nevertheless, it is interesting to observe in Fig. || that a satisfactory 
approximate size effect can be obtained by truncating the small-size asymptotic series expansion [Eq. jig)] after the 
linear term: 



a N = 




where the contants a max and L a will be discussed in what follows. Indeed, Eq. (|2(j) allows to obtain the transition 
between a horizontal asymptote characterizing the strength theory for which there is no size effect, and a decreasing 
asymptote, corresponding to a power law of exponent: —1/2 + (C — (ioc)/2 for large structure sizes. A possible 
justification is that, in the large-size first-order expansion [Eq. (Uq)], the second term in numerator has no influence 
on the power law un ~ £ _1 / 2 +(C-C!o<!)/2 f or i ar g e s i zes b u t induces a divergence at small sizes. 

One limitation of the approximate size effect is that in Eq. (p0|), the values of gm and cm'> and the values of L\ 
and L-x, are assumed to be equal and respectively characterized by a max and L a . Indeed, the estimate of these values 
would be surely different from large size data or from small size data. Thus, the approximate size effect only gives 
the shape of the size effect relation on the nominal strength but does not allow for a determination of the parameter 
values u rn ax and L a . Only the crossover length Lc and the scaling exponents (i oc and £ are univocally determined 
from the roughness analysis. However, it seems reasonnable to assume that the difference between au and cm' will 
be masked by the scatter of the experimental data. On the other hand, it seems difficult to compare the crossover 
lengths L\ and L-i (characterized by L Q in Eq. (|20|)) to the crossover length Lc using their analytical expressions, 
because the former is deduced from a mechanical approach and the latter from a roughness analysis. Nevertheless, 
from a physical point of view, both have the same meaning. For small structure sizes, i.e. L <C Lq or L oy the energy 
released by the structure is negligible, while for large sizes, i.e. L 3> Lc or L Q , the energy released is dominant. Hence, 
it seems reasonnable to assume that both crossover lengths are of the same order of magnitude (it is the case in Fig. 
H where it is assumed that L a = Lc)- 



IV. EXPERIMENT 



In the following the experimental setup of mode I fracture tests on a quasibrittle material, wood, is described. 
Two wood species have been tested: Maritime pine (Pinus pinaster Ait) and Norway spruce (Picea Abies L.). The 
average oven dry specific densities were respectively 0.55 and 0.40 and the moisture content of all specimens was 
measured between 11 and 13 %. Tests were made on modified tapered double cantilever beam specimens (TDCB) 
oriented along the longitudinal-tangential directions of wood (Fig. ||). Six sets of geometrically similar specimens 
characterized by their width L — 7.5, 11.3, 15.0, 22.5, 30.0 and 60.0 mm have been used. A straight notch is machined 
with a band saw (thickness 2 mm) and prolonged a few millimeters with a razor blade (thickness 0.2 mm). Fracture 
is obtained through uniaxial tension at a constant opening rate. The fracture surfaces are generated along the average 
longitudinal-radial plane pf wood. Load-deflection values were continuously recorded during the tests. 

The roughness analysistJ, performed from these fracture tests, has shown an anomalous roug hening [Eq. (0)] of the 
crack surfaces driven by the global and local roughness exponents, ( = 1.35 ± 0.10 and Q oc = 0.88 ± 0.07 in the case 
of pine, and, £ = 1.60 ± 0.10 and Q oc = 0.87 ± 0.07 in the case of spruce. Hence, the experimental i?-curve behaviors 
and the size effects on critical energy release rates and nominal strengths should be described respectively by Eq. (|^), 
Eq. (|) and Eq. (|o|). 
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From an equivalent linear elastic approach (which is the frame of the model described in Sec. |T| and Sec. Ill) where 
the crack lengths can be estimated from the unloading compliance of the specimens , the elastic energy release rate 
Gr are computed from the load-deflection curves for any crack length increment Aa. Two examples of the energy 
release rate evolution Gr as a function of the crack length extension Aa obtained for both wood species and for 
a specimen size L = 11.3 mm are given in Fig. ^. An i?-curve behavior is observed for both wood species., i.e. a 
pronounced evolution of the resistance to crack growth as a function of the equivalent crack length incrementu. After 
a characteristic propagation distance (Aa sa t — 23 mm for pine (a) and 40 mm for spruce (b)) which corresponds 
approximately to the ultimate load P u , the energy release rate becomes independent of the crack length increment 
Aa. The post i?-curve behavior arises at constant resistance to crack growth and corresponds to critical resistance 
Grc- In Figure || the i?-curves and post i?-curve behaviors are fitted with the Eq. (||). These fits are obtained in 
keeping three free parameters: the specific surface energy 7, the ratio A/Iq 1 ^ 100 and the scaling exponent (C — Cloc)/z. 
The expected i?-curve behavior [Eq. (||)] provides a good description of the increase of the experimental resistances to 
crack growth. For both wood species, the fitted exponents (C — Ooc)/ z are close to those measured from the roughness 
analysisE3 which are given in brackets in Fig. ^. Note that the i?-curve is more pronouced in pine than in spruce 
(i.e. the critical energy release rate Grc is greater for a smaller crack length increment Aa SQ t than in spruce). This 
trend has also been observed for the other specimen sizes L. 

On the other hand, for both.. wood species, it has been shown that the sizes L of tested specimens were greater 
than the crossover lengths Lcl3 [defined in Eq.(^) and Eq. ( po|) relative to the size effect]. Thus, according to Eq. 
([|) and for L Lc, the size effect on the critical energy release rates of both wood species is expected to evolve as 
a power law: Grc ~ L^~^ loc . The critical energy release rates obtained for both wood species are plotted in Fig. |5| 
versus the characteristic specimen sizes L. Experimental size effects are fitted by a power law Grc(L) ~ L@ whose 
the exponents f3 = 0.42 for pine (a) and (3 = 0.64 for spruce (b) are in fair agreement with the expected exponents 
C — C/oc = 0.47 ± 0.17 for pine (a) and 0.73 ± 0.17 for spruce (b). Note that the size effect on critical energies of spruce 
is stronger than the one obtained for pine. 

In the same way, in Fig. [}[ the nominal strengths <jn obtained for both wood species are plotted versus the 
characteristic specimen sizes L. According to Eq. ( |20| ) and, as previously mentioned for L Lc, the size effect 
on nominal strengths is expected to evolve as: <jn ~ i^ 1 / 2 +(C-Cioc)/2 ]\j te that latter asymptotic behaviour is 
obtained with the assumption L = Lc in Eq. (20) as previously discussed in Sec III C . As shown in Fig. ^, simple 
power law fits of the experimental nominal strengths give exponent which are in good agreement with those expected, 
i.e. — i + *'~|'° c = —0.27 ± 0.09 for pine (a) and —0.14 ± 0.09 for spruce (b). By comparison of the size effects on the 
critical energy release rates and on the nominal strengths it is appeared that, as expected intuitively, the more the 
size effect on critical energies is important the more the size effect on nominal strengths is weak. 



V. CONCLUSION 



From an energy-based analysis, a link between the scaling laws describing the fracture surfaces and the size effect 
on the nominal strength of the structures is proposed. On the basis of a Family- Vicsek scaling, which has been shown 
to reflect a purely elastic brittle fracture behaviortS, the size effect obtained is in agreement with the classical size 
effect of linear elastic fracture mechanics: .-a n ~ L^ 1 / 2 [Eq. (^)], In the case of an anomalous scaling, reflecting the 
fracture behavior of quasibrittle materialsE3, an asymptotic analysis allows to estimate the size effect relation on the 
nominal strength, and especially, the crack length increment for which the maximum load is reached. From large-size 
and small-size asymptotic series expansion of the size effect, an approximate size effect [Eq. (p0[)] of general validity 
is proposed. This relation represents a smooth transition from the case of no size effect, for small structure sizes, to a 
power law size effect which appears weaker than the size effect of LEFM. Thus, in the case of quasibrittle materials, 
the approximate size effect relation [Eq. (p0|)] is different from the classical size effect law proposed by Bazantc3 
and especially for large structure sizes where an asymptotic behavior un ~ L^/s+CC-Cioc)/ 2 j g p re( jicted instead 
of the size effect of LEFM suggested by Bazant's theory. The difference can be explained by the fact that, for an 
anomalous scaling of fracture surfaces, the critical resistance to crack growth at the maximum load evolves as a power 
law Grc ~ Z>~£ ,t>c for large structure sizes while this resistance is assumed to be constant in LEFM. Experiments 
performed on geometrically similar wood specimens of various sizes, for which anomalous roughening of crack surfaces 
has been observed previouslyEJ, show that the size effects on the critical resistances and on the nominal strengths are 
in fair agreement with the predicted asymptotic behaviors Grc ~ L^~^ loa and cat ~ 

L -i/ 2 +(C-Cio (! )/2_ on the other 

hand, if one considers a weakly anomalous roughening of the fracture surfaces {i.e. £ — > Qoc), -R-curve [Eq. (^)] and 
size effect on the critical resistance [Eq. (0)] vanish progressively and as a consequence, the size effect characteristic 
of a quasibrittle material [Eq. (pp[)] tends to the classical size effect of a purely elastic brittle material [Eq. (ph]. 
Experiments on different kinds of quasi-brittle materials are currently being performed to test our predictions. 
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FIG. 1. Size effect on the critical resistance to crack growth Grc obtained in the case of an anomalous scaling (£ = 1.3, 
Ooc = 0.8, A — 0.1 and l = 1: arbitrary values 

FIG. 2. Geometrically similar Tapered Double Cantilever Beam (TDCB) fracture specimens of different sizes 
L = 7.5, 11.3, 15.0, 22.5, 30.0 and 60.0 mm. 



FIG. 3. Approximate size effect on the nominal strength (solid curve, Eq. (|20|)) and asymptotic series expansions (dashed 
curves, Eq. @ and Eq. @) 

FIG. 4. Examples of .R-curves Gr(Ao) respectively obtained for pine (a) and spruce (b) from specimens of characteristic 
size L = 11, 3 mm. 

FIG. 5. Size effect on critical energy release rate Grc respectively for pine (a) and spruce (b). The expected slopes from the 
roughness analysisEj are ( — £; oc = 0.47 ± 0.17 for pine (a) and 0.73 ± 0.17 for spruce (b). 



FIG. |&. Size effect on nominal strength respectively for pine (a) and spruce (b). The expected slopes from the roughness 
analysis^ are - \ + c ~^'° c = -0.27 ± 0.09 for pine (a) and -0.14 ± 0.09 for spruce (b). 
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Figurc2-Morel et al.-PRB 
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